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Effect of transient pinning on stability of drops sitting on an inclined plane 
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We report on new instabilities of the quasi-static equilibrium of water drops pinned by a hy- 
drophobic inclined substrate. The contact line of a statically pinned drop exhibits three transitions 
of partial depinning: depinning of the advancing and receding parts of the contact line and depinning 
of the entire contact line leading to the drop's translational motion. We find a region of parameters 
where the classical Macdougall-Ockrent-Frenkel approach fails to estimate the critical volume of the 
statically pinned inclined drop. 



PACS numbers: ???? 

I. INTRODUCTION 

Dispense a drop on a flat substrate and then tilt it. 
Depending on the balance between gravitational and cap- 
illary pinning forces, the drop will slide down or stay at 
rest. Raindrops sticking or sliding on a vehicle windshield 
provide a familiar example of this drop stability problem, 
which is of broad practical importance. 

In structural genomics, for example, protein crystals 
are grown by dispensing protein-containing drops onto 
horizontal glass or plastic substrates. Because protein 
crystals are extremely fragile, the substrate is then in- 
verted so as to prevent any nucleating crystals from sed- 
imenting onto and adhering to it Crystals in these 

" hanging drops" instead sediment to the liquid-air inter- 
face, where they can be easily extracted without damage. 
It was shown previously [i| that pinning conditions are 
important for maintaining stability of the inclined drops 
of protein solutions with concentration suitable for crys- 
tallization. Obviously, crystal nucleation and growth - 
and thus the ultimate quality of the molecular structure 
determined by X-ray crystallography - are strongly af- 
fected by the drop shape and consequently by the drop 
stability. Drop shape variations are also a major ob- 
stacle to automated optical recognition of the protein 
drop's contents, important in high-throughput experi- 
ments. More generally, the motion of contact lines is 
related to motion of elastic manifolds in the presence of 
disorder ^ , including motion of interfaces in porous me- 
dia and depinning of flux line lattices |^ , Wigner crystals 
and charge-density waves Q . 

The pinning of an inclined drop that prevents its con- 
tinuous motion is related to the contact angle hysteresis, 
whose magnitude is usually estimated from the maxi- 
mum difference between the contact angles 9a and 9r at 
the advancing (downhill) and receding (uphill) edges of 
the contact line, as shown in Fig. ^ If this maximum 
difference A{cos9)r,a — cos9r — cos 6*0 is nonzero, then 
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FIG. 1: A drop on an inclined surface, characterized by the 
drop volume V and diameter d, the contact line perimeter p, 
the advancing and receding contact angles Oa and Or, and the 
substrate inclination angle a. 

drops of volume less than a critical VJa) may remain at 
rest at a given inclination angle a It^,^, although this is 
a necessary but not sufficient condition for drop stability. 
A simple formula 

A{cos9)r,a = AVs'ma (1) 

was obtained by Macdougall and Ockrent as a phe- 
nomenological explanation of their experiments with the 
inclined drops 9] and independently by Frenkel as a 
boundary condition of one exactly formulated problem 
of the drop stability on a tilted surface f^. Here A = 
d^^a^^ is an appropriately scaled material constant, 
a = i'^/pgY^^ is the capillary length, and p, 7, g, and d 
are the density, surface tension, gravity, and drop width, 
respectively. 

The MOF formula Eq. ^ is believed to describe the 
relation between contact angle hysteresis and the equilib- 
rium and criticality of an inclined drop 0,11, 0,0, 0,0, 
[TM [T9I I , from small inclinations and volumes at which the 
drop deforms but remains static up to the critical incli- 
nation or volume at which it begins to slide continuously. 
The validity of the MOF functional form has been veri- 
fied in a variety of experiments 8, 9, lOj JT, JX ]3e ■ 
However, surprisely it was noted that all phenomenolog- 
ical improvements have led to the same general form of 
Eq. with different length factors instead of d and 
this formula has not always correctly predicted the exact 
volumes of drops pinned at critical conditions. 
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Despite extensive study, the problem of the stabihty 
of a one-component drop on an inchned surface has yet 
to be addressed in its full richness. Here we examine the 
equilibrium and criticality of a water drop pinned on a 
hydrophobic flat glass slide. In particular, we examine a 
wide range of quasi-static tilting 0° — 90° for a moderately 
hydrophobic surface (~ 90° contact angle) having a mod- 
erate range of contact angle differences ~ 30° — 90°. We 
show that the MOF criterion is not general and its limi- 
tation is based on the existence of transient modes of dis- 
placement, whose mechanics have not yet been described. 
We have clearly resolved the partial and global instabil- 
ities of the contact line. Three coexistence curves corre- 
sponding to the partial depinning of the advancing and 
receding parts of the contact line and to depinning and 
continuous motion of the entire contact line mark those 
instabilities unambiguously. The transient displacements 
of the advancing and receding parts of the contact line 
have different scalings. We find that the difference be- 
tween the drop's maximum and minimum contact angles 
cannot reliably be used to predict the maximum pinning 
strength and the onset of drop sliding. 
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FIG. 2: Views normal to the substrate plane of the static 
contact lines for two drops, achieved as the substrate incli- 
nation is increased. The coordinates are in pixels, with the 
same image magnification for both drops, and the scale bar is 
2.6 mm long. The numbers on each contact line indicate the 
substrate inclination a in degrees. 



II. CONTACT LINE STABILITY AND 
DEPINNING FOR AN INCLINED DROP 

A drop's contact line, between its initially static pinned 
state and its steady sliding motion, may exhibit a con- 
tinuous series of intermediate states. This fact has been 
previously noted [H IH El 113 , but its effect on the sta- 
bility of an inclined drop has not been fully appreciated. 

After being dispensed on a homogeneous, flat, horizon- 
tal substrate, a drop will have a circular contact line. As 
the substrate angle a is slowly stepped upward, a drop 
of volume V eventually becomes unstable and slides con- 
tinuously down the substrate. At a smaller angle, the 
drop's contact line may become locally unstable, and un- 
dergo local displacements that change the contact line's 
shape but that do not produce continuous motion. Al- 
though these critical angles for the onset of global and 
local instability are in general different, for some exper- 
imental conditions they may be weakly distinguishable 
in measurements. In this case, the transient pinning of 
the contact line does not affect drop criticality and may 
be neglected. The drop equilibrium and criticality may 
then be described in terms of a simple energetic balance, 
as in the MOF formula, between the potential and cap- 
illary energies of the inclined drop. When the global 
and local equilibriums are well separated, the contact 
line can be displaced over the substrate while simulta- 
neously maintaining its stability against continuous slid- 
ing. In this case, because the contact line configuration 
corresponding to the global equilibrium is unknown, the 
simple energetic balance describing the drop stability cri- 
terion should be reconsidered. 

Several attempts to analytically describe the shape, 
reconfigurations and criticality of the contact line for an 



inclined drop have been reported. In 1948, using a vari- 
ational technique to analyze the drop shape for small in- 
clinations, Frenkel ^ explicitly showed for a 2D inclined 
drop that equilibrium conditions and translational drop 
instabihty lead to the MOF criteria Eq. (P) . Popova [13 
extended the variational technique to a 3D drop at small 
inclination. She analytically calculated the equilibrium 
drop shape, contact line shape, and the contact angle as 
a function of p osition along the contact line. Carre and 
Shanahan used the ideas similar to Popova's varia- 
tion of the contact angle along the contact line to calcu- 
late the pinning force, and obtained a criticalit y eq uation 
similar to the MOF criterion E g. IT l). Dussan |21| . stim- 
ulated by earlier experiments [wL |Tl|, studied critical- 
ity of 3D inclined drops with an initially elliptically-like 
contact line. Using the equations of continuous fluid dy- 
namics, she found an equation of equilibrium states for 
the inclined drop for small hysteresis, which can be used 
for obtaining the critical conditions of the inclined drops 
also. Popov [l^l used a variational analysis to exam- 
ine the equilibrium and criticality of a weakly perturbed 
hemi-spherical drop at large inclinations. His solution 
describes well only the stability of small drops inclined 
near a = 90°. 

It is important to note that in all the above studies 
0, lL6, 2^ J the contact line shape was either as- 
sumed arbitrarily or else determined from the minimum 
of some free energy. Furthermore, in all cases the cho- 
sen contact line shape was assumed to be maintained up 
to and including the critical point for the onset of insta- 
bility and drop sliding. Consequently, these approaches 
ignored the transitional behavior of the contact line prior 
to global criticality. 

In addition to these analytic attempts, two numerical 
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studies 01 nil have explored freely displacing contact 
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FIG. 3: Total displaced length of the contact line L at (a) the 
advancing and (b) the receding edges as a function of drop 
inclination angle a. The initial drop perimeter at a = 0° is 
pa, and at acrit the contact line can no longer find a static 
configuration and begins to slide continuously. The numbers 
on each curve denote the drop volumes in ^l. 

lines. Dimitrakopoulos and Higdon Esj reported tran- 
sient contact line behavior similar to that observed here. 
However, the implementation of pinning in the numeri- 
cal algorithm, which leads to the experimental-like con- 
tact line profiles, was arbitrary, and their j/-constrained 
boundary conditions need further clarification and justifi- 
cation, lliev j2J| used two phenomenological parameters 
to describe the pinning. However, the unclear connec- 
tion between those parameters and experiment and the 
absence of numerically calculated the drop deformations 
and criticality curves make comparison with the present 
results difhcult. 



III. MATERIAL AND EXPERIMENTAL 
METHODS 

Distilled and deionized water purified by a NANOp- 
ure H system (Barnstead, Boston, MA) was dispensed 
onto siliconized flat glass slides with diameter 22 mm 
(HR3-231, Hampton Research, Laguna Niguel, CA). On 
a freshly unpackaged slide, a 40 jil water drop formed 
a reproducible contact angle of 90 — 92°. To determine 
drop stability on a given slide, a drop was manually dis- 
pensed onto a horizontal slide using a 100 fil micropipette 
(Pipetman Co., France). The slide was then slowly ro- 
tated in 2 — 4° steps on a home-built rotation stage. The 
time interval between rotations was roughly one minute, 
long enough to allow transient shape relaxations to dis- 



sipate. A 640x480 pixel resolution digital camera (Cohu, 
San Diego, CA) with a telecentric 55 mm objective (Com- 
putar, Japan) was mounted on the rotation stage. Image 
recording at six frames per second began immediately 
after each stage rotation was completed and continued 
throughout the entire relaxation period. A custom im- 
age recognition program was written and implemented 
in Lab View to process each image to extract the contact 
line. Fig. |3 shows examples of contact line profiles at 
different tilt angles determined in this way. 

The apparent contact angles at the advancing and re- 
ceding contact lines were measured from the drop side 
view using an independent horizontal goniometer. Dis- 
pensed drop volumes were accurate to 0.1 — 0.5%, and tilt 
and contact angle measurements were accurate to 1 — 2°. 
Measured velocities U of average drop contact line mo- 
tion of the transient displacements were < 0.1 mm/s. 
Using water's dynamic viscosity rj ~ 0.01 g/(cm s) and 
surface tension 7 ~ 70 dyn/cm yields an upper bound 
< 10^^ for the capillary number Ca — U-q/^. Thus all 
dynamic effects during contact line rearrangement can be 
neglected. 

On a horizontal, homogeneous flat surface, a drop's 
minimum free-energy configuration has a circular contact 
line. In practice, the actual contact line shape depends 
on the initial contact conditions formed while the drop 
is dispensed. We found that the subsequent contact line 
displacements depend on the contact line's initial shape 
and on the initial contact angles along it. Consequently, 
we carefully prepared and selected drops with initially 
circular contact lines. 



IV. RESULTS 

Fig. 121 shows typical results for the contact line posi- 
tion recorded at different tilt angles, for two drops with 
volumes of 13 fjd and 90 fjd, respectively. As the substrate 
is tilted, the contact line remains pinned in its original 
circular configuration. Beyond a first critical tilt angle 
Qfa, the advancing portion of the contact line becomes lo- 
cally unstable and displaces in an attempt to find a new 
equilibrium, eventually reaching a new static configura- 
tion. Beyond a second, larger critical angle Ur, the reced- 
ing part of the contact line becomes locally unstable and 
displaces, but the drop again finds a new static configu- 
ration. Beyond a third critical angle of inclination acrit, 
the drop becomes unstable and slides continuously. The 
difference in the behavior of the advancing and receding 
contact lines implies that pinning along the contact line 
is not homogeneous. This conclusion is consistent with 
previous calculations and measurements [l^ ITsL IT^ IT^ 
of the contact angles along a contact line's circumference. 

Fig-Elexamines how the displaced length L of the static 
contact line evolves with respect to the critical parame- 
ters a and V. The total displaced length L of the contact 
line at the advancing and receding edges following a tilt 
increment to angle a can be determined by subtracting 
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images acquired at a = 0° from those at that angle a. 
For smah a such that the contact hne remains circular, 
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FIG. 4: Contact line transitions as a function of substrate 
inclination angle a and drop volume V . Symbols denote ■: 
onset of quasi-static displacements of the advancing contact 
line, □: onset of quasi-static displacements of the receding 
part, and onset of continuous motion of the whole contact 
line. 



differences in the transient displacements of the contact 
lines. Between absolute stability (zone (a)) - where the 




FIG. 5: Log-log representation of Fig. |1] The letters denote 
zones where (a), the contact line is absolutely stable against 
any inclination; (b) the contact line is stable up to a maximum 
inclination < 90°; (c) the contact line locally displaces to a 
new static configuration; (d) the contact is globally displaces 
to a new static configuration; (e) the contact line is unstable 
and moves continuously. 



the resultant images still show some displaced pixels that 
are randomly distributed over the contact line, arising 
from noise and other measurement errors. In Fig. |31 (a) 
and (b), this regime corresponds to the horizontal "zero" 
parts of the curves. At larger a (beyond aa or a^), the 
advancing or receding part of the contact line begins to 
displace (Fig. 13 insets), and the difference image shows 
a chain of connected pixels. This chain grows on further 
inclination to form a displacing front of length L. 

Fig. 13 shows the ratio of the total displaced length L 
to the unperturbed drop perimeter po (at a = 0°) versus 
Oi/oicrit 1 where acrit corresponds to onset of drop instabil- 
ity and continuous translational motion. Fig .13 vi sualizes 
unambiguously the fact noted previously |lOllllL ITM |20| 
that the advancing contact line may displace at a lower 
inclination angle than the receding contact line. Al- 
though there is considerable scatter in the data, there is 
still remarkable consistency in behavior over the factor- 
of-10 volume range examined. The advancing contact 
line begins moving at a small a /acrit < 0.2 and the 
displaced length L/po grows monotonically, reaching a 
consistent value of ~ 0.6 ± 0.1 just before Ucrit- In 
contrast, the receding contact line remains pinned un- 
til a /acrit > 0.6 — 0.8, and then steeply increases to 
L/po ^ 0.5 at acrit- Note that, because the drop be- 
comes distended, the total contact line length near acrit 
exceeds its initial length pq. 

Fig. 0] presents a subset of the data in Fig. El plotted 
in the space of critical parameters V and a. The solid 
■ and open □ squares indicate the onset of local insta- 
bility at the advancing and receding parts of the contact 
lines, respectively, and the solid triangles ▲ indicate the 
onset of continuous drop motion. A log-log scaling of 
this data is presented in Fig. [3 and clarifies the observed 



initial circular contact line is maintained at all inclina- 
tions - and continuous motion (zone (e)), the contact 
line passes through three transitions: instability of the 
advancing line at aa(V) (Fig. ^ curve ■), instability of 
the receding line at ar{V) (Fig. 0] curve □), and finally, 
instability at acritiV) leading to continuous translational 
motion of the entire contact line (curve A). Fig. [3 also 
clearly shows that there are five zones of behavior in the 
space (y, a) in which the contact line loses its stability. 
In addition to zones (a) and (e), in zone (b) the contact 
line is stable only up to a maximum inclination < 90° . In 
zone (c), between the I and □ transition curves, the con- 
tact line shows partial instability at its advancing edge. 
In zone (d) both the advancing and receding portions 
participate in quasi-static displacements, but the drop 
still remains at rest. 

It is interesting to note the functional dependencies of 
the transitional curves. The receding contact line tran- 
sition and the transition to continuous sliding have the 
same volume scaling y-o 'i's ^ However, the ad- 

vancing transition scales as V^^ '~'^ ^ V~^. This sug- 
gests that the advancing instability is controlled by dis- 
turbances that do not scale like the drop size. The length 
factor d in the MOF formula for the advancing part of 
the contact line should thus be replaced by some new 
scaling (5, which is independent of the drop size and is 
likely related to the length scale of local perturbations. 

Fig. 13 shows the measured apparent advancing and re- 
ceding contact angles versus inclination angle, for drop 
volumes ranging over a factor of 8. The measurements 
were performed for each particular drop volume at incli- 
nations below acrit, for which the contact line reached 
a static configuration after each inclination increment. 
Solid guide lines denote the stable regions ((a) and (b) 
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FIG. 6: Apparent values of the advancing contact angle 6a 
(open symbols) and receding contact angle 6r (closed sym- 
bols) vs. inclination a, for five different volumes V. The lines 
are guides to the eye, and are comprised of a solid segment 
where the contact line remains pinned in its initial a = 0° 
configuration (regions (a) and (b) in Fig. and a dashed 
segment where the contact line undergoes local displacements 
(regions (c) and (d) in Fig. 0. The contact line as a whole 
remains in static equilibrium until the end of the dashed line, 
and slides continuously beyond it (the A-line in Fig. |SJ . 



in Fig. |3J), and dashed lines indicate zones of partial in- 
stability at the advancing and receding edges ((c) and 
(d) in Fig. O. The drop contact line traverses the sta- 
ble zones (a) and (b), passes the advancing and receding 
displacement transitions and the zones of partial insta- 
bility (c) and (d), and eventually reaches the absolutely 
unstable zone (e). These data deviate significantly from 
those presented in Ref. [2^, where a hydrophobic sub- 
strate was also used. In particular, while the advancing 
angle in Fig. is nearly independent of drop volume, the 
receding angle depends strongly on drop volume. 

Fig. El raises another interesting question, touched on 
in |27l | : which values of the advancing and receding con- 
tact angles 0a and 9r for a given drop volume do we have 
to choose for an adequate description of contact line pin- 
ning and stability? According to Fig. [51 and Fig. jSl the 
advancing and receding angles often lie in different zones 
of stability in Fig. jSl In particular, the last absolutely 
stable of 9r (indicated in Fig. jSl by the points at which 
the solid curves connect to the dashed ones) do not have 
corresponding absolutely stable values of 0a ; the advanc- 
ing line has undergone quasi-static displacements that 
modify (reduce) 6a from what would be obtained if the 
advancing line had remained in its initial a = 0° posi- 
tion. In this case when the angles lie in different zones of 
stability, they cannot be simply inserted into the MOF 
formula to get an accurate measure of pinning strength 
and drop stability. Traditionally, the contact angle hys- 
teresis at criticality is obtained from the difference 9a — dr 
measured at the last inclination before the drop begins to 
slide. These correspond to the last points on the dashed 
curves in Fig.jSl and to the ▲ - transition curve in Fig.^ 
and Fig. [3 
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FIG. 7: Equilibrium and criticality of an inclined drop. Sym- 
bols indicate the evolution of drops from static pinning to 
continuous sliding, for five different volumes. The solid line 
1-2 is a fit to the MOF formula (Eq. Q), the dashed lines 
show how the data deviates from this fit and the solid curve 2- 
4 indicates instability of the contact line (the A-line in Fig.l^. 
The letters and inserts represent the same regions plotted in 
Fig. |S] and different kinds of displacement of the contact line 
respectively. For definitions of and Fg see the text. 

Fig.Qshows a direct test of the MOF equation as a pre- 
dictor of the equilibrium and stability of inclined drops. 
Using the perimeters p from images similar to those pre- 
sented in Fig. 121 and the traditional choice for the con- 
tact angles 0a and 9r to quantify the contact angle differ- 
ence, we may rescale the data of Fig. 01 and Fig. |B1 using 
= j{cos9r — cos9a)p/2TT and Fg = sin a such that 
the MOF fit of Eq. appears as a straight line. At small 
values of the scaled variables, the data are in fact linear, 
and a fit gives 7 = 72 dynes/cm, consistent with the ac- 
cepted value for water at T=22C of 72.5 dynes/cm. The 
part 1-2 of the solid line corresponds to the regions (a) 
and (b) in Fig. jsl where the contact line of pinned drops 
is in equilibrium and does not change its initial config- 
uration (see an insert corresponding to (a, b)-region in 
Fig. [3). Thus, the MOF approach is able to describe 
well deformation of a meniscus upon the drop inclination 
while the contact line remains completely pinned and the 
gravity and capillarity forces can balance each another. 
This behavior in this region of parameter space has been 
observed in many previous studies 0, 0, 0, 0, 0| . 

After the point 2 at larger Fg and F^, corresponding 
to larger inclinations or drop volumes or contact angle 
differences 9a — 9^, the equilibrium curves deviate sub- 
stantially from the linear fit: line 2-3 (dashed curves in 
Fig. Cjl. The line 2-3 corresponds to the region (c) in 
Fig. where the advancing part of the contact line is 
not stable and exhibits quasi-static displacements (see 
an insert corresponding to the region (c) in Fig. jJI). The 
points at largest Fg - on the " hooks" of the dashed curves 
in Fig. [3 - represent the last stable contact line config- 
urations, corresponding to the ▲ -line in Fig. |S1 The 
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zone between the curves 2-3 and 2-4 corresponds to the 
(d) region of Fig. \E\ According to Fig. [3 in this zone 
both the advancing and receding contact hnes are un- 
stable but only undergo quasi-static displacements (not 
sliding). These last points form the criticality curve 2-4, 
which obviously, cannot be described by the MOF for- 
mula. 



V. CONCLUSIONS 

The present experiments show that the standard pic- 
ture of instability of inclined drops is over-simplified, and 
that the behavior is more complex than assumed in previ- 
ous analytical treatments. In particular, the assumption 
that the contact line remains unchanged with inclination 
until it begins to slide does not adequately describe the 
actual depinning. 

As the substrate tilt is increased, three distinct depin- 
ning transitions - of the advancing portion of the contact 
line, the receding portion of the contact line, and of the 
contact line as a whole - are observed. These transitions 
have different effects on overall drop stability. Transient 
displacements of the advancing contact line in zone (c) 



decrease drop stability locally (B-curve in Fig.^J. But on 

further inclination, in zone (d), transient displacements 
of the receding contact line (occurring with additional ad- 
vancing displacements) involve whole contact line in re- 
configuration that stabilizes the overall pinning (A-curve 
in Fig. ^ and the "hooks" in Fig. d). We understand 
this stabilizing effect as arising a new mechanism of dis- 
sipation directed against the drop sliding, which origin is 
based on pinning-depinning of the local parts of the drop 
contact line. 

Consistent with the MOF formula, the contact angle 
difference A(cos0)r,a varies linearly with sin a up to in- 
clinations at which the receding contact line displaces. 
But the MOF formula fails to describe drop stability at 
larger angles and the critical transition to drop sliding. 
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